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We theoretically study reflection and transmission of light 
in a one-dimensional disordered phase-conjugating medium. 
Using an invariant imbedding approach a Fokker-Planck equa- 
tion for the distribution of the probe light reflectance and ex- 
pressions for the average probabilities of reflection and trans- 
mission are derived. A new crossover length scale for local- 
ization of light is found, which depends on the competition 
between phase conjugation and disorder. For weak disorder, 
our analytical results are in good agreement with numerical 
simulations. 

PACS numbers: 42.65.Hw, 78.20.Ci, 71.55.Jv, 72.15.Rn 



Over the last two decades scattering of light from 
random optical media has received a lot of attention 
PI . In passive random media many interesting multiple- 
scattering effects were discovered, such as enhanced 
backscattering of light |0], intensity correlations in re- 
flected and transmitted waves |^ and Anderson local- 
ization y. Also absorbing or amplifying random op- 
tical media have been investigated. In the latter, the 
combination of coherent amplification and confinement 
by Anderson localization leads to amplified spontaneous 
emission and laser action without using mirrors, which 
have been observed in laser dyes B and semiconductor 
powders H resp. These being linear random media, it 
is interesting to ask what happens in a nonlinear active 
random medium, such as a disordered phase-conjugating 
medium (PCM). A PCM consists of a nonlinear optical 
medium with a large third-order susceptibility x i see 
Fig. n^. The medium is pumped by two intense counter- 
propagating laser beams of frequency ujq . When a probe 
beam of frequency luq + S is incident on the material, a 
fourth beam will be generated due to the nonlinear po- 
larization of the medium. This conjugate wave has fre- 
quency ujQ — S and travels with the reversed phase in the 
opposite direction as the probe beam H. The medium 
thus acts as a "phase-conjugating mirror". Depending 
on the characteristics of the PCM, the reflected beam is 
either stronger or weaker than the incoming one, while 
the transmitted probe beam is always amplifled |^. It 
has been shown that phase conjugation also occurs in 
disordered x media |9[]. This raises several interest- 
ing questions with respect to reflection and transmission 
of light at such a disordered medium: (1) how are the 
amplifying properties of a transparent PCM affected in 
the presence of disorder? (2) What are the fundamental 
similarities and differences between a nonlinear random 



phase-conjugating medium and a linear amplifying or ab- 
sorbing random medium? (3) Is there a regime in which 
Anderson localization occurs, and what are the require- 
ments to observe this? These questions and their answers 
form the subject of this paper. 

Our starting point is the wave equation describing a 
one-dimensional (ID) disordered PCM ||l3] 
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Here /c^^ = {ujq ± S)/c and *(a;) = {£p{x),£*{x)), with 
£p{x) and £*{x) the slowly- varying amplitudes of the 
probe and conjugate electric fields respectively. The off- 
diagonal parameter 7 = 70 e"^ = -^x'"^^^i^2 is the 
pumping-induced coupling strength between the probe 
and conjugate waves in the PCM, with £1, £2 the electric 
field amplitudes of the two pump beams. The disorder 
is modeled by a randomly fluctuating part e{x) in the 
relative dielectric constant lllll. 







FIG. 1. A ID disordered phase-conjugating medium, as de- 
scribed in the text. Solid (dashed) horizontal lines denote 
probe (conjugate) beams. 

In order to calculate the reflection and transmission 
coefficients Tp, r^, t„ and tr,, we use an invariant imbed 



ding approach ||l^,|l3[. Following Ref. |l^ we obtain the 
evolution equations for the probe and conjugate waves in 
the medium 
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ko = ujq/c and /? = y^P+T^/c. Using the boundary 
conditions from Fig. ^ at x = and x = L then yields 
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In the absence of phase conjugation, for 70 = 0, equa- 
tions ^^) and ( pq ) reduce to the weU-known imbed- 
ding equations for a Unear random medium p3| , and 
Vc = tc = 0. In the absence of disorder, equations ( pq ) 
and (pq) reduce to the evolution equations for r^ and tp in 



a transparent PCM, and rp = t^ = 0. Equations (pi) sat- 
isfy the energy conservation law Rp + Tp — R^ — T^ = 
1, with Rp = |rp|2 the probe reflectance etc. They 
form the basis of all our results here. We first de- 
rive a Fokker-Planck (FP) equation for the probability 
distribution of Rp. We set rp = Rpe^^^, substitute 
this into (|5a|), subsequently into the Liouville equation 
dQ 

dL 
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where Q{Rp,Qp) is 

the density of points {Rp,Qp) in phase space, and aver- 
age over the disorder. Assuming a gaussian distribution 
for e(L), with (e(L)} = and (e(L)e(L')) = g5{L - L), 
where pointed brackets denote an average over disorder, 
yields 



dW 



Rp[l — -Rj 



\2 d^W 



{1 - 2{3 + MiL))Rp + 5Rl) 
M{L) - 2Rp)W, W{0) = 6{Rp). 



when ^0 ^ I//?- Equation (^) is of the same form as 
the equation for the probability distribution of the re- 
flectance at a linear active random medium ||lj|, with 
the important difference that in the latter case M(L) is 
an L-independent constant, proportional to the imagi- 
nary part of the dielectric constant. Using this anal- 
ogy, our phase-conjugating medium alternates between a 
linear amplifying (for M{L) < 0) and linear absorbing 
(for M{L) > 0) random medium. For M{L) = the 
well-known FP equation for a passive random medium 
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il5| is retrieved. 



(1 - Rp)^W 

Multiplying both sides of (^ by Rp and integrating by 
parts leads to a recursion relation for the moments of the 
probe reflectance, 

j^{R;) = n^R;+') 2n{n M{L)){R;) + n^R^-^). 

(7) 

For n—1 and setting (i?p) « {Rp) jl^], integration of (|^) 
yields for the average probe reflectance 

{Rp) ^[C + o?-il cos(2/3L) -I- 2aii-il sin(2/3L) - (C-h 
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with C = (4/32 + a2)(2,52 + 72) and a = l/Co- In the 
absence of phase conjugation, this reduces to {Rp) = 1 — 
e^^/So p^ a,nd in the absence of disorder {Rp) = 0, as for 
a transparent PCM. Using equations ( pa| ) and (pq), one 
can directly obtain an evolution equation for the average 
Zn^m = (ii^'i?™), which is given by 

-7jZn,m = n{m + n)Z„+i,„,. - [2?!^ + m{n + l)]2'„.„j + 
dl 

n^ Zn-i,rn + 2(m + n)M{L)Zn,m + 2mM{L), (9) 

and equivalent to (M) for m = 0. Solving (0) for the 
conjugate reflectance yields 



{Re) = -1 
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e-"^ + {Rp). (10) 



Similarly, one obtains for the probe transmittancc from 
(H) and di) PI 



{Tp) = 
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Here W = (Q), I = L/^o, M{L) = ^oB{L) and ^0"^ = 
■kgkp , the inverse localization length in the absence of 
phase conjugation. In deriving (g) we have neglected an- 
gular variations of W , the random-phase approximation 
(RPA), which applies to the situation of weak disorder 



The conjugate transmittancc is then given by {Re) — 0, 
/o\ through the conservation law {Rp) + {Tp) — (i?c) — {T^.) = 

In order to test these analytical predictions we have 
carried out numerical simulations. Using a transfer ma- 
trix method [|9|, equations (|]) are discretized on a ID 
lattice with lattice constant d, into which disorder is in- 
troduced by letting e{x) randomly fluctuate from site to 



site. Figures (0)-(0) show the probe and conjugate re- 
flectance and transmittance as a function of the length 
L of the medium for various values of the detuning 5 
and disorder. In all cases we took d — 10^* m and 
0^0 — lO^^s"^ and typical PCM parameters. 

Fig. H shows how the periodic behavior of {Re) and 
{Tp) which is characteristic of a transparent PCM be- 
comes "modulated" by an exponentially decaying enve- 
lope in the presence of weak disorder. Simultaneously, 
and with the same periodicity, some probe light is now 
reflected and some conjugate light transmitted, due to 
normal reflections in the disordered medium. When the 
amount of disorder is increased, the oscillatory behav- 
ior of the reflectances and transmittances is less and be- 
comes suppressed for large L, see Fig. g. The reflected 
probe and conjugate intensities then both saturate, with 
limL^oo(i?c) = limi,^oo(-Rp) - 1, and {Tp) and {T,,) de- 
cay to zero (localization). For a transparent PCM the 
conservation law Tp — Re — 1 applies, i.e. for each pump 
photon scattered into the forward (probe) beam in the 
medium, a photon from the other pump is scattered into 
the backward (phase-conjugate) beam. In the localiza- 
tion regime of Fig. 0, on the other hand, the conserva- 
tion law {Rp) — {Re) — 1 applies (cf. Eq. (|lO|)). Hence 
(Tp) has exchanged roles with Rp due to disorder: all 
pump photons which are absorbed into probe and conju- 
gate beams are now reflected and despite ampliflcation, 
transmitted intensities are suppressed. This suppression 
has also been found in linear amplifying random media 
1^,^. The saturation of {Re) suggests that the phase- 
conjugate reflected beam arises in the region into which 
the probe beam penetrates and that ampliflcation takes 
mostly place within a localization length of the point of 
incidence. The behavior of the transmitted intensities 
with increasing length of the medium is determined by 
two competing effects: on the one hand, enhancement 
occurs due to increased probability of multiple reflec- 
tions. On the other hand, less light is transmitted due 
to increased probability of retroreflection of the incoming 
probe light. For small L, the latter effect dominates {Tp) 
in Fig. y. As L increases, the increasing amplification of 
probe light due to multiple scattering takes over, which 
leads to exponential increase and a maximum in {Tp). For 
again larger L, most of the probe light is reflected, and 
{Tp) decreases exponentially to zero, as in a normal dis- 
ordered medium [g2|. The crossover length scale Le be- 
tween exponential increase and decrease is given by the 
solution of S^ + 7^cos2(/3L) = 27^/3^0 cos(/3L) sin(/3L), 
which for (5 ^ 70 becomes 
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dominates over phase conjugation and has almost washed 
out the oscillatory behavior of {Rp) and {Tp). 

Comparing the numerical results with the analytic ones 
from (H), (10) and ( [ll| ) we flnd good agreement (de- 
viations < 5 %) for weak disorder as in Fig. ||, and 
for stronger disorder and weak phase conjugation as in 
Fig. H. In the intermediate regime, for ^0 > 1//3 « c/70 
results differ considerably, see inset in Fig. ^. There the 
RPA and the assumption {Rp) « {Rp) are not valid and 
a different approach is needed. 




FIG. 2. Probe and conjugate reflectance and transmittance 
(averaged over 10000 realizations of the disorder) at a disor- 
dered phase-conjugating medium as a function of the length L 
(in units of the lattice spacing d) of the medium. The dashed 
curve denotes Tp in the absence of disorder. Parameters used 
are S = 10^°s"\ 70 = 4 • 10^°s"^ and e{x) G [-0.05,0.05], 
corresponding to a localization length ^0 = 1.2 m. 
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FIG. 3. Same as Fig. ti but now for larger disorder, 
e{x) £ [—0.5,0.5], corresponding to a localization length 
^0 = 12 mm. The inset compares the analytic results (H) 
for {Rp) (dashed curve) and (Ull) for {Tp) (thick solid curve) 
with the numerical ones. 



In the opposite limit of (5 S> 70 phase-conjugate reflection 
is weak (maximum value of Re — 0.16) and we retrieve 
exponential localization, see Fig. 0. Randomness now 
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FIG. 4. Same as Fig. |, but for S = 10"s-\ The thick 
dotted curves denote (Rp) and (Tp) in the absence of phase 
conjugation. 

In conclusion, we have studied reflection and trans- 
mission of light at a ID disordered phase-conjugating 
medium in the limit of small disorder, for (3^q ^ 1. 
The predicted behavior of reflectances and transmit- 
tances arising from the interplay between amplification 
and Anderson localization displays similar features as 
that in a linear disordered amplifying medium. The 
main difference is the coupling of two waves in the PCM, 
which leads to additional interference effects. In future 
work we intend to: (1) investigate the strong disorder 
regime. There the reflection of the pump beams cannot 
be neglected, and a full nonlinear analysis is required. 
(2) Study the distribution of reflection and transmission 
eigenvalues and the statistical fluctuations in reflectance 
and transmittance for a multimode 2D or 3D disordered 
phase-conjugating medium 123]. This is relevant to ex- 
periments, which mostly employ 3D PCM's pj], and in- 
teresting in the context of random lasers: in a linear 
disordered amplifying medium the average reflectance be- 
comes infinitely large with increasing amplification, upon 
approaching threshold p^. It would be interesting to 
investigate whether something similar occurs in a disor- 
dered PCM, this being a "naturally" amplifying medium 
and feasible candidate for nonlinear random lasing. 
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